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ABSTRACT 

We solve the Bethe-Salpeter equation in order to determine the spectrum of pseudoscalar and vector 
meson bound states for light as well as heavy quarks. The fermion propagators are obtained by 
solving the Schwinger-Dyson equation consistently with the Bethe-Salpeter equation, a procedure 
necessary for demonstrating the Goldstone nature of the pion in the chiral limit. Our results agree 
qualitatively and quantitatively with expectations both from current algebra for light quarks and 
from composite models for heavy quarks. 
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1. Introduction 

There has been considerable theoretical study of the strongly interacting 
bound states. However most of the work is based on non-relativistic approxima- 
tions to the bound state equations. This approach should be accurate if we consider 
the bound states which only involve heavy quarks, but is expected to fail for light 
mesons. In particular it does not take into account the fact that the light pseu- 
doscalar mesons are almost Goldstone particles and that the light quark masses are 
mostly dynamical. Furthermore in most of the studies there is no clear connection 
between the constituent quark mass which is fed into the bound state equation 
and the current algebra mass. In order to correctly describe the light mesons it is 
necessary to go to a field theoretic framework. In the present paper we study the 
quark-antiquark bound states by solving the Bethe-Salpeter (BS) equation consis- 
tently with the Schwinger-Dyson (SD) equation for the quark propagator, in such 
a way that we get zero mass pseudoscalar states in the limit of zero bare quark 
masses. To simplify the problem we work in the ladder approximation. We employ 
for our calculations theoretically and phenomenologically motivated models for the 
gluon propagator. 

The BS equation in the ladder approximation can be written as, 

S-\q + HP)X(P, q)S^(q (1 - Op) = ~i J l,x(p, k^G^k - q)-^ , (1) 

The quark propagators are determined by solving the SD equation, also in the 
ladder approximation, 

S- 1 {q)=i-m a (K c )+ l J ^ Sa (k) lv G^(k - q)-0^ , (2) 

where A c is the ultraviolet cutoff. The gluon propagator, G^ v , is modelled in terms 
of several parameters which can be fitted to experimental data. It is given in the 
Landau gauge by G tlv (k) = -(g^ - k^k u /k 2 )G(k 2 ) where, 

G(fc 2 ) = ^ , ( 1 + b l ^^±m + (27r??) V*) (*) + «(1 - u:k 2 /k 2 )e^^ 

3 [k 2 ln(x +x) V ln(x + x) J 

and x = k 2 /K 2 QCD . This model goes to the two loop form of the running coupling 
at large momentum and leads to an approximately harmonic oscillator potential in 
three dimensional configuration space, a form which is necessary to get the realistic 
spectrum of heavy mesons. The BS wave function for the pseudoscalar bound state 
can be expressed as, 

x(p,q) = 75[xo +^xi +^X2 + bM%3] ( 3 ) 

with similar decomposition for other spin and parity states. The resulting BS 
equation is simplified by expanding the wave functions, \% i n terms of Tschebyshev 
polynomials, 

Xl (V,Mi,cos^ =J2x ( r ) (q 2 ,M 2 B y^(co S e) , (4) 
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where we have set p 2 = -M%, the bound state mass squared. We keep only the 
leading order polynomial for our calculation. The contribution due to the dropped 
terms was estimated to be small. Details of the calculation for pseudoscalars are 
given in Ref. [1]. 

2. Results and discussion 

Here we describe some qualitative and quantitative features of the results for 
pseudoscalar and vector mesons. The detailed numerical results for pseudoscalar 
mesons and the parameter choices are given in Ref. [1]. The mass functions for 
different quarks are displayed in Fig. 1. The asymptotic behavior of the mass 
functions and the wave functions agrees with the one found on the basis of operator- 
product-expansion considerations. For light mesons our results satisfy the Gell- 
Mann, Oakes and Renner relation, M 2 b = [m (A c ) + m&(A c )] < qq > ch /fl h , where fh 
is the bare mass of the quark defined in Eqn. 2 and the subscript ch means that 
the quantity is computed by using m(A c ) = 0. Both m and the condensate < qq > ch 
depends on the ultraviolet cutoff A c but their product was found to be insensitive 
to A c . For heavy quarks, the relationship between M ab and [m o (0) + m 6 (0)] is roughly 
linear, in agreement with nonrelativistic limit expectations. 

We have also obtained preliminary results for vector mesons ground states. 
The overall fit was found to be better with parameters choices, a=-[32.9 MeV]~ 2 , 
ry=270 MeV with the remaining parameters same as in Ref. 1. The results for 
pseudoscalar mesons were found to be within a few % of the ones given in [1] . The 
results for some of the vector mesons are as follows: m p = 787 MeV (770 MeV), m K » 
= 978 MeV (892 MeV), m = 1250 MeV (1019 MeV), m D . = 2075 MeV (2010 MeV), 
m B , = 5325 MeV (5324 MeV), m J/4 , = 3245 MeV (3097 MeV), ™ T = 9500 MeV (9460 
MeV). The numbers in parenthesis are the experimental results. These results were 
obtained by keeping only the dominant term in the invariant decomposition of the 
vector meson wave function. More complete results will be given in a forthcoming 
publication. 

Finally we display our preliminary results for the electromagnetic and isospin 
splittings for pseudoscalar and vector mesons. Our results for the pseudoscalar 
mass splittings decrease significantly as we go from the Kaon to the B meson, in 
reasonable agreement with experiments. Quantitatively we find, K° - K + = 4.02 
MeV (input), D+ - D° = 3.95 MeV (4.77 MeV) B" - B+ = 1.0 MeV (0.1 MeV). 
For the vector meson we do not get as strong a decrease as is displayed by the 
experiments. The results are as follows: K*° - K*+ = 11 MeV (6.7 MeV), D*+ - D*° 
= 10.8 MeV (2.9 MeV) b*°-B*+ = 7.4 MeV (?) We are currently examining different 
models and calculating corrections to our result to determine if the agreement can 
be improved. 

In conclusion, we have presented a covariant treatment of qq bound states 
which is applicable for both light and heavy mesons. Qualitatively our results are 
in good agreement with current algebra results for light quarks and nonrelativistic 
limit expectations for heavy quarks. Further tests of the approach described here, 
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as well as parametrization of a s (q 2 ), will be done by obtaining a more complete 
spectrum and by computing the electromagnetic and weak form factors of these 
mesons. 
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